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Abstract 

We have theoretically designed three families of the half-metallic (HM) antiferromagnets (AFM), 
namely, La^4VOs06, La^4MoTcOg and LaAMoReOg (A = Ca, Sr, Ba), based on a systematic ah 
initio study of the ordered double perovskites \i&ABB'0% with the possible B and B' pairs from 
all the 3d, 4d and 5d transtion metal elements being considered. Electronic structure calculations 
based on first-principles density-functional theory with generalized gradient approximation (GGA) 
for more than sixty double perovskites L&CaBB'OQ have been performed using the all-electron full- 
potential linearized augmented-plane-wave method. The found HM-AFM state in these materials 
survives the full ab initio lattice constant and atomic position optimizations which were carried 
out using frozen-core full potential projector augmented wave method. It is found that the HM- 
AFM properties predicted previously in some of the double perovskites would disappear after the 
full structural optimizations. The AFM is attributed to both the superexchange mechanism and 
the generalized double exchange mechanism via the B fo g ) - O (2p 7r ) - B' fo g ) coupling and the 
latter is also believed to be the origin of the HM. Finally, in our search for the HM-AFMs, we find 
LaylCrTcOg and LavlCrReOg to be AFM insulators of an unconventional type in the sense that the 
two antiferromagnetic coupled ions consist of two different elements and that the two spin-resolved 
densities of states are no longer the same. It is hoped that our interesting predictions would 
stimulate further experimental searches for the HM-AFMs which have so far been unsuccessful. 

PACS numbers: 71.20.-b, 75.10.Lp, 75.50.Ee, 75.80,+w 
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I. INTRODUCTION 



Half-metallic (HM) ferromagnets (FM) was first discovered by de Groot, et al.\ based 
on their band structure calculations for magnetic semi-Heusler compounds NiMnSb and 
PtMnSb. Shortly afterwards, other magnetic materials such as Fe^O^- and CrO^ were 
also found to be half-metals. Half-metallic materials are characterized by the coexistence 
of metallic behaviour for one electron spin and insulating behaviour for the other. Their 
electronic density of states is completely spin polarized at the Fermi level, and the conduc- 
tivity is dominated by these metallic single-spin charge carriers. Therefore, half-metallic 
materials offer potential technological applications such as a single-spin electron source, and 
high-efficiency magnetic sensor a 4 ' 5 ' 6 . 

In this work, we search for another kind of HM materials with first-principles calcula- 
tions. In a HM material, the spin magnetic moment per unit cell is quantized, i.e., an integer 
number times Bohr magneton (fiB) - 1]5]6]7 It may occur that this integer for some HM ma- 
terials is zero. This situation has been called half-matellic antiferromagnetism . 6 i 8 i 9 i 10 Most 
properties such as full spin-polarized conduction electrons, zero spin susceptibility, and no 
Stoner continuum, of these HM antiferromagnts (AFM) are the same as those of the HM- 
FMs discussed abovefi. However, there is one important difference: a HM-AFM produces 
no macroscopic magnetic field. Therefore, a HM-AF could support 100% spin polarized 
charge transport without any net magnetization. Furthermore, since there is no symme- 
try operation (translation plus spin flip) that connects its spin-up and spin-down bands, a 
HM-AFM is qualitatively different from a conventional AFM such as bcc Cr. Because of 
their unique properties mentioned above, HM-AFMs have recently attracted considerable 
attention. Furthermore, they could be used as, e.g., a probe of the spin-polarized scanning 
tunneling microscope without perturbing the spin character of samples. The HM-AFMs are 
expected to play a vital role in the future spintronic devices which utilize the spin polariza- 
tion of the carriers. 

The first HM-AFM was proposed by van Leuken and de Groot£ on the basis of the 
Heusler compound "V^MnFegSbyln. Due to the complexity of this material, considerable 
effort has been spent to come up with other candidate materials, which may be easier to 
synthesize. In particular, based on extensive first-principles band structure calculations for 
double-perovskite-structutre oxides F^BB'Oq where B and B' are transition metal ions, 
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Pickett^ proposed the cubic double perovskite L^VMnOg to be a promising candidate for 
a HM-AFM. In this case, V and Mn have antiferromagnetically aligned magnetic moments 
and they exactly cancel each other. This finding was borned out by the assumption that 
both Mn and V ions are trivalent and more importantly, all Mn 3+ ions are in a low spin 
state of S = 1 (t\ t t\ a [)■ More recently, Park, et a/.— suggested thiospinel systems 
Mn(CrV)S4 and Fen.sCun s (VqsTii .a) S4 as possible HM-AFMs, based on their electronic 



structure studies. In Ref. |13|, the mixed-cation double perovskites LaA'VRuOg {A' = Ca, Sr, 
and Ba) were suggested to be candidates for HM-AFs. Androulakis, et a/.— have synthesized 
La2VMn06 samples which, however, has cubic, partially ordered double perovskite structure. 
Furthermore, it exhibits ferrimagnetic (FIM) behavior rather than antiferromagnetic one, 
with Mn and V being trivalent and with Mn being in high spin state. Therefore, to date, 
there has been no successful experimental realization of the HM-AFM, though it has been 
speculated in Refs. 11 3 that ordered mixed-cation double perovskites AA'BB'Oq would 



be promising candidates for the HM-AFM. 

Motivated by the theoretical and experimental works mentioned above, we have taken 
a rather thorough search for HM-AFMs among the mixed-cation double perovskites 
L&ABB'Oq. Unlike the previous theoretical effor t 11113 , we have not only performed elec- 
tronic structure calculations for a fixed crystal structure but also carried out full structural 
optimizations including lattice constants and atomic position relaxations. We believe that 
the structural optimizations would be important. In particular, though L^VMnOg was 
predicted to be a HM-AFM in Ref. Q, it was found to be ferrimagnetic in a partly ordered 
double perovskite structure^. We have explored a variety of pairs BB' with B and B' from 
3d (V, Cr, Co, Ni), Ad (Mo, Tc, Ru, Rh) and 5d (W, Re, Os, Ir) transition metals. Fortu- 
nately, we find double perovskites LaAVOsOg and L&AMoYOq (Y = Tc, Re) to be robust 
HM-AFMs. 



II. THEORY AND COMPUTATIONAL DETAILS 

Our ab initio theoretical search was based on first-principles calculations within density- 
functional theory (DFT)^ with the local density approximation (LDA) plus the general- 
ized gradient corrections (GGA) 1 - . We used the highly accurate all-electron full-potential 
linearized augmented plane wave (FLAPW) method 1718 as implemented in the WIEN2k 
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package^ 9 -, to calculate the total energy, electronic band structure and magnetic properties 
for the materials with fixed structural parameters. The FLAPW method makes no shape 
approximations to the electron density or potential and retains high variational freedom in 
all regions, and hence it is well suited to open crystal structures with low site-symmetries 
such as those considered here. The wave function, charge density, and potential were ex- 
panded in term of the spherical harmonics inside the muffin-tin spheres. The cut off angular 
momentum (L max ) of 10 used for the wave function and of 6 used for the charge density 
and potential are sufficient for accurate total-energy calculations 19 . The wave function out- 
side the muffin-tin spheres was expanded in terms of the augmented plane waves. A large 
number of augmented plane waves (about 115 per atom) (i.e., R mt K max = 6) were included 
in the present calculations. The set of basis functions was supplemented with local orbitals 
for additional flexibility in representing the valence V 3d states, Ru Ad states, Os 5d states 
and La 4/ states, as well as the semicore Ca (or Sr, Ba) 3s, 3p, La 5s, 5p, O 2s states. The 
muffin-tin sphere radii used are 2.5 a.u. for La and Ca (or Sr, Ba), and 2.0 a.u. for V, Ru, 
Os, Mo, Cr, Tc and Re, 1.4 a.u. for O. The improved tetrahedron method is used for the 
Brillouin-zone integration^. We used 140 ^-points in the irreducible Brillouin zone wedge, 
which correspond to 1500 Appoints in the first Brillouin-zone. 

For the structural optimization calculations to determine theoretical lattice constants 
and atomic positions, we used the faster frozen-core full-potential projector augmented wave 
method (PAW)2i as implemented in the VASP packag A cutoff energy of 450 eV for plane 
waves is used. A 8x8x6 Monkhorst-Pack fc-point grid in the Brillouin zone was used, which 
correspond to 30 /c-points in the irreducible Brillouin-zone wedge. Their atomic positions 
and lattice constants were fully relaxed by a conjugate gradient technique. Theoretical 
equilibrium structures were obtained when the forces acting on all the atoms and the stresses 
were less than 0.05 eV/A and 1.5 kBar, respectively. 

III. CRYSTAL STRUCTURE AND THEORETICAL SEARCH STRATEGY 

We consider L&ABB'O e (A = Ca, Sr, Ba) in an ordered double perovskite structure 
[space group P4/nmm (No. 129)], as shown schematically in Fig. 1. This structure can be 
regarded as a combination of cubic Lal?03 and AB'0 3 perovskites in a superlattice having 
B and B' layers stacked along the [111] direction, as in La2CrFe06- Transition metal ions 
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FIG. 1: (Color online, the file size of this figure is too big to be included) The crystal structure of 
an ordered double perovskite (P4/nmm). 

B and B' can couple to each other either ferromagnetically or antiferromagnetically in the 
[111] directions. Each B (B ; ) ion has six B' (B) neighbors. LaS03 and AB'0 3 perovskites 
can also be stacked along the [001] direction, resulting in a double perovskite structure with 
the space group P4mm (No. 99), as has been considered by Park, et alM. In this case, each 
B ion has four B ion and two B' ion neighbors, while each B' ion has four B' ion and two 
B ion neighbors. We have performed electronic structure calculations for LaylVRuOg and 
LaAVOsOe in both P4/nmm and P4mm structures. We find, however, that LaAVRuOe and 
LaAVOsOg in the [001] stacked structure are nonmagnetic and ferrimagnetic, respectively. 
Furthermore, the calculated total energy for the [001] stacked structure is about 0.25~0.35 
eV/f.u. higher than that of the [111] stacked structure. This shows that the [111] stacked 
structure is more stable and thus we will not consider the [001] stacked structure further. 

The [111] stacked structure (P4/nmm) contains twenty atoms per unit cell, i.e., two 
chemical formula units (f.u.). It has a tetragonal symmetry with lattice constants a and 
c. The atomic positions of two La and two A atoms are (1/4,3/4,0), (3/4,1/4,0) and 
(1/4,3/4,1/2), (3/4,1/4,1/2), respectively. The atomic positions of two B and two B' 
atoms are (1/4, 1/4, 1/4 + 5), (3/4, 3/4, 3/4 - 5) and (3/4, 3/4, 1/4 + 5'), (1/4, 1/4, 3/4 - 5'), 
respectively. The O atoms are divided into three different types. The atomic positions of 
two type one and two type two O atoms are Oi: (1/4, 1/4, zx), (3/4,3/4,1 — z\) and O2: 
(1/4, 1/4, 1/2 — 22), (3/4, 3/4, 1/2 + Z2), respectively, while the atomic positions of eight type 
three O atoms are 3 : (1/2 — x, 1/2 — x, 1/4 — z 3 ), (x, x, 1/4 — z 3 ), (x, 1/2 — x, 1/4 — 2 3 ), 
{l/2-x,x, I/4-23), (l/2 + x,l-x, 3/4 + ^), (l-x,l/2 + x, 3/4 + 23), (l-x,l-x, 3/4 + 23), 
(1/2 + x, 1/2 + x, 3/4 + 2 3 ). Therefore, in this structure, apart from lattice constants a and 
c, there are six more parameters to optimize, namely, 5, 5', x, zi, 22, 23. 

For the chemical formula L&ABB'Oq, B and B' can be an any pair from the twenty-four 
transition metal elements, as shown in Fig. 2, and A is one of the three alkali earth elements 
Ca, Sr, and Ba. Therefore, although we have settled down on the ordered double perovskite 
structure, there are still 3xC| 4 = 828 possible compounds to consider. This is an impossible 
task. Nevertheless, our preliminary calculations and also the existing literature tell us that 
Ti, Cu, Zr, Nb, Pd, Ag, Hf, Ta, Pt and Au in the double perovskite structure are usually 
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FIG. 2: (Color online) Part of the Periodic Table showing 3d, Ad and 5c? transition metals. Bold 
italic characters denote the twelve elements that have been selected as the BB' pairs for Ij&ABB'Oq 
in this work. 



not magnetic and thus, they can be excluded from our further consideration. Furthermore, 
we note that in the double perovskite structure, Mn and Fe are almost always in a high 
spin state (see, e.g., Refs. Consequantly, it is not possible to find a high spin state 

element from the remaining twelve elements to match Mn and Fe, and hence Mn and Fe are 
excluded from further consideration too. Therefore, in our first stage of search for HM-AFs, 
we consider all the possible BB' pairs from only the twelve transition metal elements as 
denoted by bold italic characters in Fig. 2. Furthermore, we initially confine ourselves to 
A = Ca. Therefore, in our first stage of search, we have performed self-consistent electronic 
structure calculations only for sixty-six (d, 2 = 66)) La.Ca.BB'0 6 compounds in the ordered 
doub le p erovskite structure with fixed lattice constants a = 5.492 A, c = \[2a taken from 
Ref. [LJ. When we find any LaCa-B-B'Og to be a HM-AFM, we replace Ca with Sr and Ba 
and perform further self-consistent electronic structure calculations. 

Fortunately, we find from our electronic structure calculations for sixty-six L&C&BB'O^ 
compounds, that six B and B' pairs, namely, VRu, VOs, MoTc, MoRe, NiTc, NiRe, would 
give rise to HM-AF L&C&BB'Oq compounds. We note that previously, Park, et alM the- 
oretically predict LaAVRuOg to be HM-AFMs. However, we performed both volume and 
atomic position optimization calculations for these systems and found, unfortunately, that 



the HM-AFM properties predicted in Ref. [13|, disappear. This perhaps explains why the 
attempt to synthesize the HM-AFM LaAVRuOe by Liu, et al., 24 was unsuccessful. There- 
fore, in our second stage of search, we carried out full lattice constant and atomic position 
relaxation calculations for all the possible HM-AFM L&ABB'Oq compounds found in the 
first stage of search. We find fortunately that among the possible HM-AFM LaABB'OQ 
compounds, LaAVOsOg, LaAMoTcOg and LaAMoReOg stay as the stable HM-AFMs. 
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FIG. 3: Total and orbital-decomposed density of states of LaSrVRuO-6 in the ideal P4/nmm 
structure (no atomic position relaxation, see text and Table I). 



TABLE I: Calculated physical properties of La^VRuOg in the ideal P4/nmm structure (S = 8' = 
X = Z\ = z 2 = z 3 = 0). 



A = 


Ca 




Sr 




Ba 




magnetic states 


FM/NM 


AF 


FM/NM 


AF 


FM/NM 


AF 


a (A) 


5.510 


5.517 


5.545 


5.553 


5.604 


5.612 


c/a 


V2 




V2 




V2 




v (A 3 /f.u.) 


118.31 


118.72 


120.57 


121.05 


124.44 


125.00 


m v (a*b) 


0.039 


-0.949 


0.017 


-0.967 


0.006 


-0.992 


mR u Ob) 


0.074 


0.730 


0.028 


0.738 


0.016 


0.744 


m t Ofi/f.u.) 


0.165 


0.0 


0.065 


0.0 


0.033 


0.0 


N(E F ) T 


2.563 


0.0 


2.575 


0.0 


2.358 


0.0 


(states/eV/f.u.) [ 


2.491 


4.064 


2.640 


4.205 


2.460 


4.414 


gap (eV) T 




0.74 




0.79 




0.79 


AE AF ~ FM (eV/f.u.) 


-0.102 




-0.103 




-0.090 
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TABLE II: Calculated physical properties of La^VOsOg in the ideal P4/nmm structure (S = 5' = 
x = z\ = z 2 = Z3 = 0). 



A = 




Ca 




Sr 




Ba 




magnetic states 




FM/NM 


AF 


FM/NM 


AF 


FM/NM 


AF 


a (A) 




5.522 


5 54fi 


5 557 


5 577 


5 607 

tj . \J\J I 


5 6S5 

•J . U O '.J 


c/a, 




V 








V - 1 




v ( A 3 /f u ) 




119.05 


120.61 


121.36 


122.67 


124.62 


126.53 


777,1/ |//R 1 




0.0 


-1.254 


0.0 


-1.274 


0.0 


-1.305 






0.0 


0.890 


0.0 


0.899 


0.0 


0.910 


m t (/xs/f.u.) 




0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


N(E F ) 


T 


3.420 


0.0 


3.462 


0.0 


2.527 


0.0 


(states/eV/f.u.) 


I 


3.420 


3.798 


3.462 


3.883 


2.527 


4.023 


gap (eV) 


T 




0.82 




0.87 




0.98 


AE AF ~ FM {eV/f.u. 


) 


-0.231 




-0.236 




-0.237 





IV. RESULTS AND DISCUSSION 

A. Half-metallic antiferromagnets from the initial search 

In the first round of search for the HM-AFMs, we have performed ab initio electronic 
structure calculations for sixty-six LaCai^-B'Og compounds in the ideal P4/nmm structure 
(i.e., 5 = 8' = x = Z\ = z 2 = z 3 = 0) with the lattice constants a = 5.492 A, c = y/2a taken 
from Ref. as mentioned before. Remarkably, we found LaCaVFOe, LaCaMoy'06 and 
LaCaNiF'0 6 (Y = Ru, OS; Y' = Tc, Re) to be HM-AFMs. We then replaced Ca in these 
compounds with either Sr or Ba and performed further electronic structure calculations. We 
found again that the resultant compounds are HM-AFMs. This result could be expected 
because A-site atoms in ABO3 perovskite are known to behave like a carrier reservoir and a 
volume conserver. The interaction between an A atom and its neighboring atoms is usually 
very weak such that many ABO3S with different As, have similar electronic properties. 

In the ionic picture, the atoms in the ordered double perovskites have the nominal valence 
states as La 3+ A 2+ (B B') 7+ . Therefore, in, e.g., LaAVRuOg, the transition metal atoms B 



S 



and B' can have valence configurations of V 3+ (3c? 2 ) and Ru 4+ (4c? 4 ), and the antiferromagnetic 
coupling of the high spin V 3+ (t 2g f, S=l) and the low spin Ru 4+ (t 2g ] t\ g I, S=l) states 
would give rise to a zero total magnetic moment, i.e., a AFM state, Alternatively, the 
transition metal atoms B and B' can have valence configurations of V 4+ (3<i 1 ) and Ru 3+ (4<i 5 ), 
and the antiferromagnetic coupling of the spin V 4+ {t\ |, S=l/2) and the low spin Ru 3+ 
(tl g | tl g I, S=l/2) states. Of course, in the real materials, the above simple ionic model 
would be modified because of hybridization between V 3d and O 2p orbitals and also between 
Ru 4c? and O 2p orbitals. The calculated total and orbital-decomposed density of states 
(DOS) of LaSrVRuOe in the AFM state are shown in Fig. 3. Clearly, LaSrVRuOe is half- 
metallic with a band gap of about 0.8 eV on the spin-up channel. The electronic structure 
consists of the O 2p dominant lower valence band between 2.5 and 8.0 eV below the Fermi 
level (Ep), the spin-up Ru t 2g dominant upper valence band between 0.3 and 2.0 eV below 
Ep and spin-down Ru t 2g -V t 2g hybridized conduction band between 2.0 eV below Ep and 
~1.0 eV above Ep as well as the spin-up V t 2g dominant, Ru e g dominant and V e g dominant 
upper conduction bands about 0.8 eV above Ep (see Fig. 3). The calculated DOS spectra 
are very similar to the ones reported in Ref. [lj. Our calculated occupation numbers are 
0.72 e (spin-up) and 1.68 e (spin-down) at the V-site, and 2.72 e (spin-up) and 1.98 e (spin- 
down) on the Ru-site. Therefore, the calculated charge configurations are V 2 6+ (3d 2A ) and 
Ru 3,3+ (4g? 4 " 7 ) , suggesting that the valence configurations of V 3+ /Ru 4+ and V 4+ /Ru 3+ are 
perhaps nearly degenerate and hence highly mixed, as pointed out before 1 ^. This notion is 
further supported by the fact that the calculated local spin magnetic moments of V and Ru 
are not integer numbers of 2 or 1 but less than 1.0 (see Table I). 

Similarly, we can understand the formation of the AFM state in the other HM-AFM com- 
pouds in the ionic picture. As for LaAVRuOe, in LaAVOsOe, the transition metal atoms B 
and B' can have valence configurations of V 3+ (3cP) and Os 4+ (5g? 4 ), and the antiferromagnetic 
coupling of the high spin V 3+ (t\ f, S=l) and the low spin Os 4+ {t\ g f t\ g |, S=l) states 
would give rise to a zero total magnetic moment. In LaAMoTc(Re)06, the transition metal 
atoms B and B' may have valence configurations of Mo 3+ (4<i 3 ) and Tc 4+ (4g? 3 ) [Re 4+ (5c? 3 )] 
with the antiferromagnetic coupling of the low spin Mo 3+ (t\ g j t\ |, S=l/2) and the low 
spin Tc(Re) 4+ (t\ ] t\ |, S=l/2) states. In Lay4NiTc(Re)0 6 , the transition metal atoms 
B and B' may have valence configurations of Ni 2+ (3d 8 ) and Tc 5+ (4<i 2 ) [Re 5+ (5c/ 2 )] with the 
antiferromagnetic coupling of the high spin Ni 2+ [t\ g | t\ g \ e 2 |, S=l) and the high spin 
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TABLE III: Calculated physical properties of La^lMoTcOg in the ideal P4/nmm structure (5 = 
5' = x = zi = z 2 = z 3 = 0). 



A = 




Ca 




Sr 




Ba 




magnetic states 




FM/NM 


AF 


FM/NM 


AF 


FM/NM 


AF 


a (A) 




5.629 


5.640 


5.657 


5.667 


5.699 


5.717 


c/a 




V2 




V2 




V2 




v ( A 3 /f u ) 




126.13 


126.87 


127.99 


128.66 


130.91 


132.16 


TYlMr, (ur) 
1 ' v ivi o y^iD } 




0.0 


-1.007 


0.0 


-1.030 


0.0 


-1.064 






0.0 


1.112 


0.0 


1.133 


0.0 


1.164 


m t (n B /i.u.) 




0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


N(E F ) 


T 


3.168 


0.0 


3.225 


0.0 


3.321 


0.0 


(states/eV/f.u.) 


I 


3.168 


3.160 


3.225 


3.184 


3.321 


3.271 


gap (eV) 


T 




0.44 




0.49 




0.60 


AE AF ~ FM {eV/f.u. 


) 


-0.146 




-0.159 




-0.171 





Tc(Re) 5+ (tl g 1 1% |, S=l) states. 

To see whether the HM-AFM is the stable state in these six compounds, we performed 
the FM electronic structure calculations. Unfortunately, we find that the total energy of the 
HM-AFM state in LaCaNiF'0 6 (Y' = Tc, Re) is substantially higher than that of the FM 
state, by 0.19 and 0.18 eV/formula unit (f.u.), respectively. This suggests that the HM-AFM 
state in these two systems is not stable and hence we will not consider them further. We 
then performed volume relaxation calculations for L&AVYOq and LaAMoF'C>6 (A = Ca, 
Sr, Ba; Y = Ru, Os; Y' = Tc, Re) starting with both the FM and AFM states. We find all 
these systems to remain as a stable HM-AFM after the volume relaxation. The calculated 
electronic and magnetic properties of these four families of the double perovskites are listed 
in Tables I- IV, respectively. 

Tables I-IV show clearly that the total energy of the HM-AFM state in these four fam- 
ilies of the double perovskites is significantly lower than that of the corresponding FM or 
nonmagnetic (NM) state. The energy difference is larger than 0.1 eV/f.u. and can be as 
large as ~ 0.24 eV/f.u., e.g., in LaBaVOs0 6 (Table II). The insulating gap in the spin-up 
band structure is rather large and can be up to 1.0 eV as in LaBaVOs0 6 (Table II). It is 
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TABLE IV: Calculated physical properties of La^lMoReOg in the ideal P4/nmm structure (5 = 
5' = x = z\ = z 2 = z 3 = 0). 



A = 




Ca 




Sr 




Ba 




magnetic states 




FM/NM 


AF 


FM/NM 


AF 


FM/NM 


AF 


a (A) 




5.640 


5.656 


5.669 


5.687 


5.739 


5.739 


c/a 




V2 




V2 




V2 




v ( A 3 /f u ) 




126.89 


127.93 


128.85 


130.06 


133.67 


133.64 


TYlMn (ur) 

' ' V 1V1 O Xf^D } 




0.0 


-1.138 


0.0 


-1.175 


0.0 


-1.210 


m Re (fj, B ) 




0.0 


1.100 


0.0 


1.133 


0.0 


1.161 


m t (/ie/f.u.) 




0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


N(E F ) 


T 


2.925 


0.0 


3.003 


0.0 


3.175 


0.0 


(states/eV/f.u.) 


I 


2.925 


2.190 


3.003 


1.853 


3.175 


1.423 


gap (eV) 


T 




0.16 




0.27 




0.35 


AE AF ~ FM {eV/f.u. 


) 


-0.148 




-0.172 




-0.214 





also interesting to note that no strong FM solution can be stabilized in all the four families 
of the compounds, and hence these compounds have either a weak FM or NM metastable 
state (Tables I- IV). Nonetheless, the local magnetic moments in the HM-AFM state are not 
large either, all being in the order of 1.0 /ie/atom. 

B. Effects of structural optimization 

In our second stage of search, we carried out full structural optimization (i.e., both 
lattice constant and atomic position relaxations) calculations for all the possible HM-AFM 
L&ABB'Og compounds described in the previous subsession. Interestingly, we find that 
LaAVOs0 6 , LaAMoTc0 6 and LaAMoRe0 6 (A = Ca, Sr, and Ba) remain to be the stable 
HM-AFMs, although LaAVRu0 6 become nonmagnetic after the full structural optimization. 
The equilibrium lattice constants and atomic positions as well as the stability of the HM- 
AFM of these three families of the double perovskites are listed in Table V. Clearly, the 
HM-AFM state is stable over the FM or NM state in these compounds, and the size of the 
energy difference between the HM-AFM and FM/NM states is similar to the one before the 
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full structural optimization (see Tables I-IV). The full structural optimization further lower 
the total energy of these compounds rather significantly. The most pronounced decrease 
in the total energy occurs in the Ba-based double perovskites. For example, the decrease 
is about 0.2~0.3 eV/f.u. for LaBaMoF0 6 (Y = Tc, Re) and about 0.3~0.4 eV/f.u. for 
LaBaVOsCV This may be attribued to the fact that the ionic radius of Ba is significantly 
larger than that of Ca. This notion is further supported by the fact that the equilibrium 
unit cell volume increases as one replaces Ca with Sr and then with Ba (see Table V). The 
decrease in the total energy become smaller as one moves from Ba through Sr to Ca. The 
decrease in the total energy is about 0.1~0.15 eV/f.u. for LaSrMoyOe and about 0.05~0.1 
eV/f.u. for LaCaMoF0 6 (Y = Tc, Re). 

The shape (i.e., the c/a ratio) of the unit cell for the ordered double perovskites in the 
HM-AFM state does not change much after the full structural optimization (Table V). In 
particular, the deviation of the c/a for L&C&BB'Qq from the ideal value of \/2 is less than 
1 %. This is because La and Ca have rather similar atomic radii of 3.92 and 4.12 a.u.— , 
respectively. The deviation increases very slightly as Ca is replaced by Sr and Ba, and 
the largest deviation of about 1.3 % occurs in LaBaVOsOe- Sr and Ba have larger atomic 
radii of 4.49 and 4.65, respectively. However, the internal atomic displacements after the 
structural optimization are rather pronounced. For example, the contraction (expansion) of 
the in-plane B-0 (B'-O) bondlength of the BOq (B'Oq) octohedra which is described by the 
parameter x, can be as large as 2.7 % (or 0.038 A) in LaCaMoTc0 6 . The in-plane B-O-B' 
bond bending which is described by the parameters Z3, S and 5', is also rather significant. 
In particular, the B-O-B 1 bond bending can be as large as 8° (i.e., the angle of the B- 
O-B' bond is about 172°) in LaBaVOsOg. A similar amount of the bond bending occurs 
in LaBaMoTcOe and LaBaMoReOe- This may be attributed to the fact that the atomic 
radius of Ba is considerably larger than that of La. This bond bending becomes smaller when 
Ba is replaced by Sr and Ca. Another significant atomic displacement is the contraction 
(expansion) of the apical B-0 (B'-O) bond of the -B06 (B'Oq) octohedra which is described 
by the parameters Z\, Z2, 5 and 5' '. For example, the contraction of the V-O2 bondlength is 
about 4.3 %. Interestingly, this displacement of the apical O atoms is asymmetrical. Indeed, 
the V-Oi bond is elongated slightly, i.e., expands by 0.21 %. 

Obviously, the full structural optimization has the largest effects on LaAVRuOe in which 
both the HM-AF and FM states disappear and the systems become nonmagnetic after the 
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FIG. 4: Total and orbital-decomposed density of states of LaSrVRuO-6 in the theoretical deter- 
mined P4/nmm structure (full structural optimization, see text). 

structural optimization, as mentioned before. The lowering of the total energy due to the 
structural optimization is the largest among the four families of the HM-AFM candidates 
found in our first stage of search, varying from 0.4 to 0.8 eV/f.u. Likewise, the contraction of 
the V0 6 octahedra and the expansion of the Ru0 6 octahedra during the structural optimiza- 
tion are also the largest (see Tables V and VI). For example, the contraction (expansion) 
of the V-0 (Ru-O) bond on the xy plane is about 3.0 % and that of the apical V-0 2 bond 
is as large as 6.4 %. Furthermore, the B-O-B' bond bending is around 7° in all the three 
LaAVRuOe compounds. It is therefore conceivable that the weak AFM and FM state of 
LaAVRu0 6 found in the ideal P4/nmm structure could not sustain such large structural 
distortions and both of them become nonmagnetic. 

Why the magnetic states disappear in LaAVRuOe after the full structural optimization 
can perhaps be best understood by examining their electronic structures. Fig. 4 shows 
the total and orbital-decomposed densities of states of LaSrVRu0 6 in the fully optimized 
structure. They are similar to the ones for LaSrVRuOe in the ideal P4/nmm structure 
displayed in Fig. 3, except that there are not spin-polarized anymore. The electronic 
structure consists of the O 2p dominant lower valence band between 2.7 and 7.7 eV below 
E F) the Ru t 2g dominant upper valence band between 0.0 and 1.9 eV below E F and V t 2g 
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dominant lower conduction band between 0.0 eV and ~1.0 eV above Ep as well as the V 
t2 9 /Ru e g dominant upper conduction bands about 0.5 eV above Ep (see Fig. 4). The 
most striking feature is the pseudo-gap at Ep and LaSrVRuOg is a semimetal with a low 
DOS of 0.3 states/eV/spin/f.u. at Ep. It is clear that because of the pseudo-gap at Ep, 
LaylVRuOg are now stable against magnetic instabilities and therefore the systems stay in 
the nonmagnetic state . 26 ' 27 In the nonmagnetic band structures of L&AVOsOq, LaAMoTcOe 
and LaAMoReOe, there is no such pseudo-gap at Ep and hence they retain the HM-AFM 
state after the full structural optimization. 

C. Electronic structure and magnetic properties 

The calculated electronic and magnetic properties of robust HM-AFM L&ABB'Oq are 
listed in Table VII. Clearly, all three families of LaAVOsOg, LaAMoTcOg and LaAMoReOg 
(A = Ca, Sr, Ba) are HM-AFMs with an insulating gap in the spin-up channel. The spin- 
down density of states at Ep is rather large in LaAVOsOg and LslAMoTcOq, although it 
is rather small in LaylMoReOg. The insulating gap in the spin-up channel is also rather 
large (0.5~1.0 eV) in LaAVOsOe and LaAMoTcOg but relatively small (within 0.5 eV) in 
LaylMoReOg. The calculated local spin magnetic moments in all the double perovskites are 
not large and in the order of 1.0 /x^/atom. This seems to be consistent with the simple ionic 
picture described in Subsec. A in the case of LaAMoTc(Re)06 but at varience with it in the 
case of LaAVOsOg. The deviation from the simple ionic model is more apparent in terms 
of the calculated occupation numbers of the transition metal d orbitals. For example, in 
LaSrVOsOg, the calculated occupation numbers are 0.72 e (spin-up) and 1.72 e (spin-down) 
at the V-site, and 2.40 e (spin- up) and 1.69 e (spin-down) on the Os-site, giving rise to 
a charge configuration of V 2 ' 6+ (3c? 2 ' 4 ) and Os 3 ' 9+ (4c? 4 ' 1 ). In LaSrMoTcOg, the calculated 
occupation numbers are 0.74 e (spin-up) and 1.82 e (spin-down) at the Mo-site, and 2.40 e 
(spin-up) and 1.23 e (spin-down) on the Tc-site, resulting in a charge configuration of Mo 3 4+ 
(3d 2 6 ) and Tc 3,4+ (4c? 3 - 6 ). In LaSrMoRe0 6 , the calculated occupation numbers are 0.68 e 
(spin-up) and 1.94 e (spin-down) at the Mo-site, and 2.14 e (spin-up) and 0.95 e (spin-down) 
on the Re-site, resulting in a charge configuration of Mo 3,4+ (3d 2,6 ) and Re 3,9+ (5d 3A ). 

The total and orbital-decomposed DOS spectra of robust HM-AFM LaSrVOsOg, 
LaSrMoTcOg and LaSrMoReOg are shown in Figs. 5, 6 and 7, respectively. As mentioned 
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TABLE V: Theoretical structural parameters and stability of L&ABB'Oq in the HM-AFM as well 
as FM and NM states. The lattice constant a is in the unit of A. 



(a) La^VOsOg 


A = 


Ca Sr 


Ba 


magnetic states 


FM/NM AF FM/NM AF 


FM/NM AF 


a (c/a) 

5(5') (nr 2 ) 

x( Zl ) (10- 2 ) 
z 2 (z 3 ) (10- 2 ) 

V(A 3 ) 


5.538 (1.401) 5.541 (1.405) 5.585(1.395) 5.580(1.402) 5.652(1.385) 5.647(1.396) 
0.76 (0.39) 0.71 (0.22) 0.46 (0.33) 0.39 (0.14) -0.09 (-0.16) -0.03 (-0.13) 
0.60 (0.42) 0.44 (0.46) 0.60 (0.35) 0.42 (0.34) 0.56 (0.22) 0.38 (0.22) 
0.79 (0.59) 0.52 (0.62) 0.98 (1.14) 0.70 (1.18) 1.34 (2.07) 0.90 (1.88) 
119.06 119.50 121.53 121.83 125.04 125.67 


AE AF-FM ( e v/f. u .) 


-0.104 -0.102 


-0.099 


(b) LaAMoTc0 6 


A = 


Ca Sr 


Ba 


magnetic states 


FM/NM AF FM/NM AF 


FM/NM AF 


a (c/a) 
5(5') (10- 2 ) 
x( Zl ) (HT 2 ) 

Z2(ZZ) (KT 2 ) 

V(A 3 ) 


5.643 (1.408) 5.656(1.408) 5.679(1.404) 5.689(1.406) 5.740(1.392) 5.747 (1.399) 
0.49 (0.08) 0.47 (0.07) 0.14 (-0.01) 0.15 (-0.11) -0.34 (0.06) -0.19 (-0.25) 
0.17 (0.26) 0.67 (0.16) 0.20 (0.15) 0.02 (0.07) 0.24 (-0.10) 0.02 (-0.02) 
0.07 (0.58) -0.10 (0.03) 0.26 (1.17) 0.01 (1.36) 0.62 (1.93) 0.16 (1.98) 
126.50 127.35 128.54 129.41 131.64 132.77 


ae af-fm ( e v/f.u.) 


-0.135 -0.149 


-0.163 


(c) LaAMoRe0 6 


A = 


Ca Sr 


Ba 


magnetic states 


FM/NM AF FM/NM AF 


FM/NM AF 


a(c/a) 
5(5') (10- 2 ) 
x(zi) (10- 2 ) 

Z2(Z 3 ) (10- 2 ) 

v (A 3 ) 


5.651(1.409) 5.671(1.408) 5.687(1.408) 5.703(1.409) 5.749(1.393) 5.756 (1.404) 
0.40 (0.17) 0.38 (0.14) 0.20 (0.03) 0.09 (-0.07) -0.14 (-0.18) -0.14 (-0.28) 
-0.08 (-0.03) -0.18 (-0.07) -0.08 (-0.05) -0.18 (-0.11) -0.06 (-0.09) -0.21 (-0.13) 
-0.13 (0.30) -0.28 (0.57) -0.10 (0.86) -0.25 (1.09) -0.04 (1.80) -0.27 (1.94) 
127.15 128.45 129.45 130.68 132.33 133.81 


ae af-fm ( e v/f.u.) 


-0.182 -0.229 


-0.271 
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TABLE VI: Theoretical structural parameters of LaAVRuOg. 



A = 


Ca 


Sr 


Ba 


magnetic state 


NM 


NM 


NM 


a (A) 


5.517 


5.559 


5.629 


c/a 


1.409 


1.403 


1.393 


6(6') (10- 2 ) 


0.90 (0.40) 


0.36 (0.36) 


-0.29 (0.23) 


x( Zl ) (10- 2 ) 


0.78 (0.52) 


0.78 (0.38) 


0.78 (0.24) 


Z2(Z 3 ) (10- 2 ) 


0.91 (0.57) 


1.25 (1.19) 


1.70 (2.06) 


V (A 3 /f-u.) 


118.26 


120.46 


124.26 



before, when Sr is replaced by either Ca or Ba, the electronic structures of the resultant 
double perovskites are very similar to that of the corresponding Sr-based double perovskite, 
and therefore they are not shown here. Comparison of Fig. 5 with Fig. 3 shows that the 
electronic structure of LaSrVOs0 6 is very similar to that of LaSrVRu0 6 in the hypothetical 
HM-AFM state, with Os 5d playing the role of Ru Ad. Therefore, as for LaSrVRu0 6 in 
the AFM state, the electronic structure of LaSrVOsOg consists of the O 2p dominant lower 
valence band between 3.0 and 8.4 eV below Ep, the spin-up Os t 2g dominant upper valence 
band between 0.2 and 2.0 eV below E F and spin-down Os t 2g -\ t 2g hybridized conduction 
band between 2.0 eV below E F and ~1.0 eV above E F as well as the spin-up V t 2g dominant, 
Os e g dominant and V e g dominant upper conduction bands about 0.8 eV above Ep (see Fig. 
5). LaSrVOs0 6 is half-metallic with a band gap of about 0.8 eV on the spin-up channel. 
Fig. 5 shows that the lower valence band contains rather discernable contributions from 
transition metal B and B' d orbitals especially Os bd in the lower part of the band, because 
of the hybridization between O p and B (B') d orbitals. This rather strong O p - B (B r ) d 
hybridization pushes the V (Os) e g dominant band above the V (Os) t 2g dominant band, as 
shown in Fig. 5, resulting in a crystal field splitting of the t 2g and e g bands of about 2.0 eV 
for the V0 6 octahedron and of a slightly larger value for the Os0 6 octahedron. This B d - 
O p - B' d (T-bonding is the well known superexchange coupling between the B and B' ions 
and gives rise to the AFM exchange interaction between the B and B' ions. However, in the 
conventional superexchange- type AFM insulators such as NiO, the majority (minority) spin 
d orbitals of each transition metal ion are separated by a large superexchange energy gap 
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from the minority (majority) spin d orbitals of its nearest neighbor transition metal ions. In 
contrast, the majority spin V t 2g and minority spin Os t 2g orbitals in LaSrVOsOe are closely 
located in energy (Fig. 5). This gives rise to a strong hybridization between the majority 
spin V t 2g and minority spin Os t 2g orbitals via O p^ orbitals which results in the broad B 
and B' t 2g hybridized spin down conduction band in the energy range of -2.0~1.0 eV, and 
further strengthens the AFM exchange coupling between the V and Os ions. Clearly, this 
coupling between the V and Os t 2g orbitals via O p^ orbitals would allow the spin-down 
conduction electrons to hop freely from a V ion to the neighboring Os ions and back to the 
V ion, resulting in a lower kinetic energy. This situation is similar to the double exchange 
mechanism— of the metallic ferromagnetism in the colossal magneto-resistive manganites, 
and hence may be called the generalized double exchange mechanism. The direct on-site 
exchange interaction would then lower the majority spin Os t 2g band below Ep and push 
the minority spin V t 2g band above Ep, thereby creating an insulating gap in the spin up 
channel (see Fig. 5). As a result, the band gap may be called an antiferromagnetic coupling 
gap 7 . Therefore, the antiferromagnetism in LaSrVOsOg has the dominant contributions 
from both the superexchange mechanism and the generalized double exchange mechanism 
(see Refs. 29 i 30 for the detailed analysis on these two mechanisms in the doped ferromagnetic 
manganites). The origin of the half-metallicity is essentially the same as that of the HM 
ferrimagnetic double perovskites such as Sr 2 FeMo0 6 , namely, the B (t 2g ) - O (2p n ) - B (t 2g ) 
hybridization, as illustrated in Refs. QJ3 It could then be infered from this discussion that 
the antiferromagnetic transition temperature T/v in the HM-AFMs would be in the same 
order of magnitude as that in S^FeMoOe (T c = 419 K)&, i.e., being above room temperature. 

The electronic structures of LaSrMoTcOg and LaSrMoReOg are similar to that of 
LaSrVOsOe, as can be seen from Figs. 5-7. In particular, the electronic structure of 
LaSrMoTcOg is nearly identical to that of LaSrVOsOg. The discernable difference in the 
electronic structure between LaSrMoTcOe and LaSrVOsOg is the larger Mo 4d contribution 
to the O 2p dominant lower valence band, indicating a stronger Mo Ad - O 2p bonding 
in LaSrMoTcOe- This may be expected because the Mo 4c? orbitals are significantly more 
extended than that of V 3d. There are also some discernable differences in the electronic 
structure between LaSrMoReOg and LaSrVOsOg. The most pronounced difference is that 
there is a small dip or pseudop gap at Ep in the spin-down channel. Therefore, the spin- 
down band structure of LaSrMoReOe may be regarded as a semimetal with a small DOS 
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FIG. 5: Total and orbital-decomposed density of states of LaSrVOsOe in the theoretical determined 
P4/nmm structure (full structural optimization, see text). 

of 0.6 states/eV/spin/f.u. at Ep. The spin-down DOS at Ep furthers decreases when Sr 
is replaced by Ba (see Table VII). Further differences among the electronic strucrtures of 
LaSrVOsOe, LaSrMoTcOe and LaSrMoReOe are minor and they include slight differences 
in the energy position and bandwidth of various bands. 

D. Unconventional antiferromagnetic insulators 

In our first round of search for the HM-AFMs, we also found LaACrTcOe and LaACrReOe 
be antiferromagnets. However, LaACrTcOg and LaACrReOg are insulators. Nonetheless, 
we may regard these materials as unconventional AFM insulators because they have several 
subtle differences in the electronic structure from archetypical AFM insulators such as MnO 
and NiO in the cubic rocksalt structure. The calculated electronic and magnetic properties 
of these materials are listed in Table VIII. The total and orbital-decomposed DOS spectra 
of LaSrCrTcOg are displayed in Fig. 8, as an example. These results are obtained from 
selfconsistent electronic structure calculations for the fully theoretically optimized crystal 
structures. Compared with the HM-AFMs, these compounds have a stronger spin splitting 
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TABLE VII: Calculated magnetic and electronic properties of LavlXYOg in theoretically deter- 
mined P4/nmm structure using GGA. 

(a) La^VOsQg 



A= Ca Sr Ba 



magnetic states 




FM/NM 


AF 


FM/NM 


AF 


FM/NM 


AF 
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(b) LaAMo 


TcO-6 








A = 




Ca 




Sr 




Ba 




magnetic states 




FM/NM 


AF 


FM/NM 


AF 


FM/NM 


AF 


mMo {hb) 




0.0 


-1.029 


0.0 


-1.088 


0.169 


-1.128 






0.0 


1.125 


0.0 


1.176 


0.120 


1.213 


m t {n B /i.\i.) 




0.0 


0.0 


0.0 


0.0 


0.463 


0.0 


N(E F ) 


T 


3.192 


0.0 


3.206 


0.0 


3.067 


0.0 


(states/eV/f.u.) 


I 


3.192 


3.210 


3.206 


3.106 


2.572 


3.003 


gap (eV) 


T 




0.52 




0.68 




0.79 


(c) LaAMoRe0 6 


A = 




Ca 




Sr 




Ba 




magnetic states 




FM/NM 


AF 


FM/NM 


AF 


FM/NM 


AF 


m mo (Ms) 




0.0 


-1.226 


0.0 


-1.264 


0.0 


-1.301 


m Re (fj, B ) 




0.0 


1.168 


0.0 


1.199 


0.0 


1.230 


m t (/ze/f.u.) 




0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


N(E F ) 


T 


2.947 


0.0 


3.086 


0.0 


3.333 


0.0 


(states/eV/f.u.) 


I 


2.947 


1.134 


3.086 


0.594 


3.333 


0.175 


gap (eV) 


T 




0.25 




0.35 




0.49 
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FIG. 6: Total and orbital-decomposed density of states of LaSrMoTcG-6 in the theoretical deter 
mined P4/nmm structure (full structural optimization, see text). 
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FIG. 7: Total and orbital-decomposed density of states of LaSrMoReOe in the theoretical deter 
mined P4/nmm structure (full structural optimization, see text). 
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TABLE VIII: Calculated structural parameters, magnetic and electronic properties of antiferro- 



magnetic insulators LaACrYOg 


in theoretically determined P4/ 


nmm structures. 






(a) La^CrTc0 6 






A = 


Ca 


Sr 


Ba 


a (A) 


5.534 


5.573 


5.638 


c/a 


1.4123 


1.4117 


1.4068 


V (A d /f.u.) 


119.66 


122.20 


126.08 


mcr (mb) 


-2.241 


-2.264 


-2.295 


m Tc (mb) 


1.683 


1.696 


1.710 


m 4 (//s/f.u.) 


0.0 


0.0 


0.0 


gap (eV) T 


1.252 


1.224 


1.197 


1 


0.653 


0.680 


0.762 


(b) La,4CrRe0 6 


A = 


Ca 


Sr 


Ba 


a (A) 


5.551 


5.583 


5.650 


c/a 


1.4117 


1.4111 


1.4070 


v (A 3 /f.u.) 


120.72 


122.79 


126.89 


"T-Cr (MB) 


-2.217 


-2.237 


-2.271 


m Tc (Mb) 


1.468 


1.484 


1.503 


m t (// B /f.u.) 


0.0 


0.0 


0.0 


gap (eV) T 


0.63 


0.76 


0.85 


I 


0.98 


1.09 


1.12 



of the B(B') t 2g band and a larger local magnetic moment on the B and B' sites (see 
Fig. 8 and Table VIII). For instance, the spin splitting of the Cr t2 9 band in LaSrCrTcOe 
is about 2.5 eV, significantly larger than that of the V t 2g band of less than 1.5 eV in 
LaSrVOs0 6 . The local magnetic moment on the Cr site in LaACrTc(Re)0 6 is twice as 
large as that on the V site in LaAVOs0 6 . Furthermore, the spin-down Cr t 2g band in 
LaACrTc(Re)06 is aligned with the spin- up Tc(Re) t 2g band, rather than with the spin- 
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FIG. 8: Total and orbital-decomposed density of states of LaSrCrTcC>6 in the theoretical deter- 
mined P4/nmm structure. 

down Tc(Re) t<i g band as is the case in the HM-AFM compounds discussed in the previous 
Subsec. Consequently, in LaATc(Re)06, the generalized double exchange mechanism is 
ineffective and the antiferromagnetic coupling is due solely to the traditional superexchange 
mechanism. This gives rise to the AFM insulating behavior in these compounds. 

In the conventional AFM insulators, the two antiferromagnetic coupled ions are made of 
the same atomic species, e.g., Ni in NiO. The two spin-decomposed total DOS spectra are 
identical and the insulating gaps for spin up and spin down channels are the same. The 
novel features of the present AFM insulators include that the two spin-decomposed total 
DOS spectra are no longer the same and the spin up insulating gap may have a size different 
from that of the spin down insulating gap, as shown in Table VIII and Fig. 8, and that 
the antiferromagnetic coupled ions consist of two different atomic species, e.g., Cr and Tc 
in LaACrTc0 6 . 
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V. CONCLUSIONS 



In conclusion, in order to search for the fascinating half-metallic antiferromagnetic ma- 
terials, we have carried out a systematic ab initio study of the ordered double perovskites 
L&ABB'Oq with the possible B and B' pairs from all the 3d, 4d and 5d transtion metal ele- 
ments being considered. Electronic structure calculations based on first-principles density- 
functional theory with GGA for more than sixty double perovskites L&Ca.BB'0 G have been 
performed using the all-electron FLAPW method. As a result, we find three families of 
the HM-AFMs, namely, LaAVOs0 6 , LaAMoTc0 6 and LaAMoRe0 6 (A = Ca, Sr, Ba). 
The found HM-AFM state in these materials subsequently survives the full ab initio lattice 
constant and atomic position optimizations which were carried out using frozen-core full 
potential PAW method. We attribute the AFM to both the superexchange mechanism and 
the generalized double exchange mechanism via the B (i 2s ) - O (2p w ) - B' (t 2g ) coupling and 
also believe the latter to be the origin of the HM. We also find that the HM-AFM prop- 
erties predicted previously in some of the double perovskites would disappear after the full 
structural optimizations, thereby providing an explanation why the previous experimental 
attempts to synthesize the HM-AFM double perovskites did not succeed. Finally, in our 
search for the HM-AFMs, we find LaACrTcOe and LaACrReOe to be AFM insulators of 
an unconventional type in the sense that the two antiferromagnetic coupled ions consist of 
two different elements and that the two spin-resolved densities of states are no longer the 
same. Stimulated by previous theoretical predictions of the HM-AFM in L^MnVOg^ and 
also in LaAVRuOe 13 , experimental effort to synthesize the HM-AFMs has been reported 
by Androulakis, et a/,—, and attempted by Liu, et alj^, respectively, which were, however, 
unsuccessful. It is hoped that our predictions of robust HM-AFMs would encourage further 
experimental searches for the HM-AFMs. 
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